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00 ' Abstract 

; ■ We investigate the consequences of long range rapidity correlations 

^^ ' in the Glasma. Particles produced locally in the transverse plane are 

^^ , correlated by approximately boost invariant flux tubes of longitudinal 

color electric and magnetic fields that are formed when two sheets of 
Colored Glass Condensate pass through one another, each acquiring a 
modified color charge density in the collision. We argue that such long 
range rapidity correlations persist during the evolution of the Quark Gluon 
Plasma formed later in the collision. When combined with transverse fiow, 
these correlations reproduce many of the features of the recently observed 
ridge events in heavy ion collisions at RHIC. 
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1 Introduction 



Among the more striking features in cosmology are the large scale fluctuations 
seen in the cosmological microwave background (CMB) and in the mass den- 
sity perturbations that result in galaxies. Such large scale fluctuations are at 
first sight hard to understand because they extend over much larger distance 
scales than could be set up by interactions of the thermal medium produced 



after the big bang. It is now believed that these large scale fluctuations origi- 
nate in small quantum fluctuations present during the inflationary epoch. Dur- 
ing the rapid expansion of the universe in this epoch, these quantum fluctua- 
tions were stretched to size scales much larger than those that were causally 
connected in the post-inflationary era when the universe was expanding in a 
state close to thermal equilibrium. Therefore such super horizon scale fluctu- 
ations cannot be much affected by the sub-horizon scale processes allowable 
in the post-inflationary thermal universe. This explains why CMB measure- 
ments provide extremely valuable information about the inflationary epoch of 
the universe, despite the fact that the CMB radiation was produced long after 
(icMB '^ 4 • 10^ years) the primordial fluctuations that are responsible for its 
features (iinfiation ^ lO^^'^ seconds). 

There is a concrete analog of such super-horizon fluctuations in the matter 
produced in high energy hadronic collisions such as heavy ion collisions at RHIC, 
as illustrated in fig. 1. In this figure, we represent the "event horizons" as seen 
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Figure 1: The red and green cones are the location of the events in causal 
relationship with the particles A and B respectively. Their intersection is the 
location in space-time of the events that may correlate the particles A and B. 



from the last rescattering of two particles A and B on the freeze-out surface. 
These are the red and green cones pointing to the past. Any event that has a 
causal influence on the particles Aor B must take place inside the corresponding 
event horizon. Any event that induces a correlation between the particles A and 
B must lie in the overlap of their event horizons. Therefore, if the particles A 
and B have rapidities y^ and j/^ , the processes that caused their correlations 
must have occurred before the time^ 
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^ We assume here that a particle detected with momentum rapidity y originates from a point 
of space-time rapidity r; ss j/ on the freeze-out surface. This is a consequence of the boost 
invariance of the coUision (at high energy), and of the fact that the local thermal motion 
spreads the rapidities by at most one unit in rapidity. 



Therefore, we see that long range rapidity correlations can only be created at 
early times-shortly after the collision or even in the wavefunctions of the incom- 
ing projectiles, that form sheets of Color Glass Condensate [l]-[3] at high ener- 
gies. In a high energy collision of these Color Glass Condensates, an interacting 
and evolving system of high intensity color electric and color magnetic fields is 
produced [4]- [12]. This collection of primordial fields is the Glasma [8,9], and 
initially it is composed of only rapidity independent longitudinal color electric 
and magnetic fields. These fields generate topological Chern-Simons charge [13]. 
Correlations associated with particle production from these fields span large dis- 
tances in rapidity. In contrast, due to the longitudinal expansion of the matter 
produced in RHIC collisions, thermal effects can only affect particle correlations 
on a distance scale of approximately one unit of rapidity. 

Because the longitudinal fields are approximately rapidity invariant, there 
are long range correlations built into the initial conditions-these inevitably have 
their origin in the quantum mechanical wavefunctions of the hadrons. Of course, 
the long range correlations are only approximately rapidity invariant. A proper 
treatment of quantum fluctuations in the hadron wavefunctions suggests that 
the characteristic distance scales- beyond which one has significant variations 
in the correlation- is of order l/as{Qs)- Here Qg is the saturation momen- 
tum of partons in the nuclear wavefunctions; it grows rapidly with both the 
energy and the nuclear size. The strong interaction strength as{Qs) is therefore 
correspondingly weak for nuclear collisions or very high energy hadron-hadron 
collisions resulting in a characteristic distance scale of several units in rapidity 
for the rapidity correlations. 

The existence of long range rapidity correlations in high energy hadronic 
collisions has been measured in ISR experiments [15], and is intrinsic to string 
models of high energy collisions such as the Lund model and the Dual Parton 
Model [16]-[17]. The essential new feature of the Glasma is that the fields are 
localized in the transverse scale over distances (of order \/Qs) that are smaller 
than the nucleon size. Because the dynamics of these fields in collisions at high 
energies or for large nuclei occurs at small transverse distances, this dynamics 
can be described by a weak coupling expansion. A further point of departure is 
that in addition to the longitudinal color electric fields envisioned int he Dual 
Parton model and the Lund model, there is a longitudinal magnetic field of equal 
intensity. It is this combination of electric and magnetic fields that generates a 
finite topological charge density. 

Such topologically distinct field configurations are of great interest in sev- 
eral areas of theoretical physics. For example, the analog of these fields in 
electroweak theory may be responsible for generating the baryon asymmetry of 
the universe [18]-[19]. In QCD, they may be the source of masses of hadrons [20]. 

Recent experimental studies have shown that there are long range rapidity 
correlations at RHIC [21]. In the STAR experiment, forward-backward cor- 
relations in the total multiplicity were studied as a function of centrality in 
Au-Au collisions. A strong correlation was found, which increases significantly 
with greater collision centrality. This correlation is stronger than expected from 
Monte-Carlo models of particle production. This effect can be understood as 



arising from the interference between the classical contribution of rapidity in- 
dependent Glasma fields, and the first order quantum correction to this result. 
The latter is short ranged in rapidity [22]. 

More specifically, the forward backward correlation as a function of the ra- 
pidities at which the multiplicity is measured is 

' dN cIn '' 



dy'-, dyL / , / 

where dN/dyi and dN/dy2 are measured in a single event, y is typically cho- 
sen to be the midpoint between j/i and 2/2 and the brackets denote the average 
over events. The denominator normalizes the expression such that C(j/,y) = 1. 
If classical fields dominated this correlation function, then C{yi,y2) = 1 since 
the initial Glasma fields are boost invariant. However, as mentioned, there is a 
short range contribution to C^g (2/1,^2) arising from quantum corrections which 
is of order a. Thus Cpg(2/i, j/2) is typically less than unity because for long 
range rapidity correlations the numerator of this expression only contains the 
long range contribution while the denominator, corresponding to zero separa- 
tion in rapidity, is a sum of both short range and long range contributions. 
For events in nuclear collisions with increasing centrality, the coupling constant 
gets progressively smaller, thereby leading to a diminished contribution of the 
short range correlation relative to the long range classical correlation. Thus 
Cpg(j/i,y2) approaches unity with increasing centrality of the collision. Such 
a forward-backward correlations, while consistent with the Glasma hypothesis, 
is perhaps not as direct a verification of the longitudinal electric and magnetic 
fields as might be desired. We will argue in this paper that so called "ridge" 
events discovered by STAR may provide a more direct confirmation of this fun- 
damental property of the Glasma. 

These striking "ridge" events were revealed in studies of the near side spec- 
trum of correlated pairs of hadrons by the STAR collaboration [23,24]. The 
spectrum of correlated pairs on the near side of the detector (defined by an ac- 
companying unquenched jet spectrum) extends across the entire detector accep- 
tance in pseudo-rapidity of order A?; ~ 2 units but is strongly coUimatcd for az- 
imuthal angles A(/). Preliminary analyses of measurements by the PHENIX [25] 
and PHOBOS [26] collaborations appear to corroborate the STAR results. In 
the latter case, with a high momentum trigger, the ridge is observed to span 
the even wider PHOBOS acceptance in pseudo-rapidity of A?7 ^ 6 units. 

A plot of the ridge in the rj — (p plane is shown in fig. 2(a). The quantity 
^pI \fPref plotted as a function of rapidity and azimuthal angle, is the density 
of particles correlated with a particle emitted at zero rapidity [43] . All particles 
with p±^ > 150 McV are included. The quantity Ap is the difference in densities 
between single events and mixed events. The quantity p^e/ comes from mixed 
events. The results are corrected for the effects of azimuthally asymmetric flow. 
As shown in fig. 2(b), an important feature of the ridge is that its height is 
strongly dependent on centrality. As the centrality of the heavy ion collision 
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Figure 2: (a) Top Figure: The ridge as seen in the measurement of two par- 
ticle correlations with minimal cut on particle momenta. Mixed events have 
been subtracted, as have been the effects of azimuthally asymmetric flow. The 
centrality bin here is 19-28% b) Bottom Figure: The height of the ridge as a 
function of the number of binary collisions per participant. Both figures are 
preliminary STAR figures from Rcf. [43] 



is increased, there is a rapid transition to the regime of long range rapidity 
correlations; there is an equally distinct, if less dramatic, coUimation of the 
width in A0 with increasing centrality, as shown in fig. 3. 
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Figure 3: The width in rapidity and azimuthal angle of the ridge as a function 
of the number of binary collisions per participant. 



These features of the ridge are not seen in proton-proton or deuteron-gold 
collisions and appear to be unique to nucleus-nucleus collisions. As mentioned 
previously, long range rapidity correlations are not unique to nucleus-nucleus 
collisions and have been seen in proton-proton collisions as far back as the 
Intersecting Storage Ring (ISR) experiments at CERN [15]. However, in that 
case, the correlations are not collimated in azimuthal angle and therefore do not 
have the striking ridge like structure observed at RHIC. 

In this paper, we argue that the ridge is formed as a consequence of both long 
range rapidity correlations that are generic in hadronic and nuclear collisions at 
high energies plus the radial flow of the hot partonic matter that is specific to 
high energy nuclear collisions. Here we deal exclusively with the ridge as seen 
in the total multiplicity of associated particles (with a minimum p± cutoff as 
low as 150 MeV). This has the advantage for us that the multiplicity density 
of particles associated with the ridge should be approximately conserved. A 
discussion of the ridge for high transverse momentum particles, while possible 
in our formalism, is complicated due to the interaction of such high momenta 
particles with the media, where there is energy loss. 

In our simple picture, the longitudinal electric and magnetic fields of the 
Glasma form flux tubes that emit a radiation spectrum isotropic in the relative 
azimuthal angle between the particle pairs. While of finite amplitude, such a 
distribution would appear featureless on a plot of the spectrum in the Arj-Acf) 
plane. The coUimation in A0 is a consequence of strong final state effects in the 
medium produced in nuclear collisions. For the particular effect we describe, 



we will see that it is a consequence of strong radial flow in the medium as also 
previously suggested in the literature [27,29,28]. Besides the strong medium 
efi^ect on the A(j> distribution, increasing the centrality of nuclear collisions also 
enhances long range rapidity correlations; these are proportional to l/al{Ql). 
Because the saturation scale is Q^ ex L, where L is thickness of the projectiles 
along the beam axis at a given impact parameter^, the average Q'^ grows with 
increasing centrality, leading to stronger long range rapidity correlations. 

Our computation of the correlated two particle distribution function for 
P-L J 9-L > Qs gives a result of the form'^ 

dN2 \ / dN \ / dN \ 

dypd^p^dyqd^q^ / \ dy^d^p^ / \ dy^d^q^ / 

1 / dN \ / dN \ 



Si^Ql \ dypd^p^ I \ dyqd'^q 

The quantity k is a constant which we shall compute explicitly. The factor of 
1/S±Q'^ has the simple physical interpretation of the area of the flux tube from 
which the particles are emitted divided by the overall area of the system"*. The 
Glasma flux tubes are illustrated in flg. 4. 



_Ap_ 
dypd'^p± I \ dy^d'-^ 



The quantity '^ is equivalent to the left hand side of eq. (3) divided by 
^ ' ^ and multiplied by the multiplicity per unit rapidity, 

'^^^ ^ S.Ql. (4) 



dy / as{Q 

where k' « 1/13.5 for an SU{i) gauge theory [45]. We therefore obtain 

A/9 K„ 



yJPref OLs 



(5) 



where we shall show that K^ is a constant of order unity. 

This relationship is basically a consequence of dimensionality-the correla- 
tions are due to a classical effect and there is only one dimensional scale which 
characterizes the Glasma. This formula is not corrected for transverse flow, 
which modifies this result, as we shall demonstrate, by introducing a depen- 
dence of the result on the azimuthal angle between the pairs. 

As mentioned, our mechanism for the ridge has features in common with the 
work of Voloshin [27] and Shuryak [28]. However, as we shall discuss, there are 
important qualitative and quantitative differences with this approach. We also 



^See Ref. [32] and references therein. 

•^Our definition of tiie differential two particle correlation is such that 

d^p^dypd^q^dyg-—- — — = A^-^ . 

dypd^p^dyqd^q^ 

■'Note that the right hand side is proportional to the area since the total multiplicity per 
unit rapidity scales as the area. 




Figure 4: Glasma flux tubes. The transverse size of the flux tubes is of order 

1/Qs. 

note that there are several other models of the ridge [30]-we will not attempt to 
discuss these here. A key difference between the dominant particle production 
mechanism we will describe here and other mechanisms described in the litera- 
ture is that our correlation is not formed from a a daughter particle splitting off 
a parent particle. Instead, the dominant QCD contribution to the long range 
rapidity correlation comes from particles that are produced independently along 
the length of flux tube of color electric and magnetic fields, localized in a region 
of size l/Qs in the transverse plane of the colliding nuclei. 



2 Computing the Glasma 2-particle correlation 

We now turn to a quantitative analysis of two particle correlations in the Glasma 
and their role in forming the near side ridge in A- A collisions. The variance of 
the two particle multiplicity distribution, for two particles with momenta p and 
q, is defined as 



C(p,q) 



dNo 



dyp(Pp±_dyq(Pq^ 



dN 



dypd^p^^ I \ dy^d'^q 



dN 



(6) 



where the brackets ( • • • ) denote an average over events. Also, see footnote 3. 
The contributions to C{p, q) have distinct origins. If we compute emission in 
a fixed configuration of the color sources, there are connected and disconnected 
pieces. These are shown in fig. 5. The top diagram is disconnected and we 
might naively think it would contribute only to the uncorrelated second term 
of the above equation. This is not true because, when we average over the 
color sources, there are contributions which involve contractions of the sources 
between the two superficially disconnected diagrams. It is these contractions 
which dominate the computation of C(p, q) because they arise from a classical 
contribution and are therefore leading order in powers of as . The bottom dia- 
gram involves an interference between a connected and a disconnected diagram 



and is suppressed by a power of the coupling^. In fact, a more sophisticated 
renormahzation group treatment [34] shows that by computing only the classical 
contribution (top diagram of fig. 5) to C{p, q), and by averaging it with evolved 
distributions of sources with rapidity, one is automatically including both the 
leading logarithmic terms of the bottom diagram as well as the leading logs 
from NLO correction to the top diagram. Therefore, in the rest of this paper, 
we consider only the top diagram of fig. 5. 
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Figure 5: Top Figure: A classical diagram which yields a non-vanishing two 
particle correlation after averaging over the color sources. Bottom Figure: A 
contribution to the correlation function associated with a quantum correction 
to the classical field. 

The classical contribution in C{p, q) can be computed analytically^ for Qg <C 
P-Lj Q-L- For p±, q± < Qg, the computation is non-perturbative and must be 
performed numerically. Techniques developed previously [5-7] to compute single 
inclusive gluon can be extended to this case. Nevertheless, it is instructive to 
compute C{p, q) in the large momentum region for several reasons. It will 
demonstrate that the effect is genuine, and we believe key features will persist 
at lower pair momenta even if numerical coefhcicnts may be different in the two 
regimes. Also, as the ridge likely persists for high momentum hadron pairs [35], 
our computation may have direct application to that case modulo the energy 
loss effects mentioned previously. 



^In the figure shown, we have kept only those diagrams which do not vanish under pairwise 
contraction of sources. There may be additional terms which involve contraction of three 
sources, but these are suppressed in simple models of the source color charge distribution and 
we do not therefore include them here. 

^It can also be computed analytically for either pj_ or qj_ 2> Qs but we shall not discuss 
that case here. 



The correlated two particle inclusive distribution can be expressed as 

C(P,9) = ^^ E {{\MTyip,qr)-{\MUp)\'){\Mi,{qr)) , (7) 



a.a' : 



A. A' 



where the classical contribution to the amplitude for the production of a pair 
of gluons with momenta p and q is 

MZ'iP,q) = e^^{p)eUq)pVA^-^p)A-'-^-'{q), 

Mlip) = e^^{p)p'A^^'^{p). (8) 

Here the e's are the polarization vectors of the gluons and a, a' are the color 
indices of the gauge fields. The average (• • • ) in eq. (7) is an average over the 
color configurations of the two nuclei; this average will be discussed further 
shortly. 

At large transverse momenta, the classical gauge fields produced in the nu- 
clear colhsion can be expressed exphcitly [4,36,37] as 

pA (p)- ^Uc ^ J ^^^^,L {p,k^) klip^-k^r ■ 

Here fate are the SU(3) structure constants, L^ is the well known^ Lipatov 
vertex [38] and pi , p2 are respectively the Fourier transforms of the color charge 
densities in the two nuclei [3]. Taking the modulus squared of eq. (8), summing 
over the polarizations and color indices, the first term in eq. (7) can be expressed 
as 

'"*'^'^'' " 64(27r)6 y-f'''"'-f^'bsfabcfa'bc} J H (-2^)2^2^ 

L^{ p,kl±)L^'{p,k2±) 

T^{p,q;{k.^) . (10) 



(p_L -fcl-L)^(P± -k21_f 

L,,{q,k:ij_)L"{q,kij_) 



The scalar product of two Lipatov vectors is 

L^{p,k^)L^{pJ^) = -\ [5^i5'"^ + e^i^i^n] k^{p^-k^)H'Ap^-l^r , (H) 
P± 

and we denote 



~.b , 



T^{p,q;{k,^) ^ip*l{k2i_)p*l{k,^)p\{k,^)p\{k 



3±J 



X P*2(P± - fc2±)p*2(g± - fe4±)P2(P± - fcl-L)/52(g± - ^3^) ) -(12) 



^Thc components of this four vector are given explicitly by L+(p, fc^) = ^ 

L-(p,fcx) = '^^'"V"'"'^ . L"(p,fc±) = -2fcY. 
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The average in eq. (12) corresponds to 



(O) 
In the MV model [1,2], 



[Dp,Dp2]W[pi]W[p2]0[pi,p2] 



W[p] = exp - / d^x_i_ 



p''{x^)p^{x^] 



(13) 



(14) 



where p can be either pi or p2. The color charge squared per unit area /i^, 
besides the nuclear radius R, is the only dimensionful scale in the problem-as 
we will discuss later, the saturation scale Qs can be expressed simply in terms 
of this scale. We will consider this Gaussian model in the rest of this paper®. 
For these Gaussian correlations, in momentum space, 



{p*\k^)p\k'^)) = {2^fpl S^'5{k^ fcl) . 



(15) 
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Figure 6: Trivial color correlation. This type of connection between the sources 
leads to a non correlated contribution to the 2-gluon spectrum, that cancels in 
the difference in eq. (6). 

Examining the structure of T in eq. (12), one observes that one of the nine 
possible quadratic combinations of the pi's and p2S is a disconnected piece, 
represented in the fig. 6, whose expression is nothing but 



(|A(p)p)(|^(9)P) 



(16) 



which is identical to the product of single inclusive distributions. It exactly 
cancels the disconnected contribution to pair production-the second term in 
eq. (7). Therefore only eight terms contribute to the correlated distributions of 
pairs. Of these, as we shall see, four terms give identical leading contributions 
to Ci(p, q) for large p, q^ Qs- Two of these terms, as shown in fig. 7(a), have a 



*In the simplest treatment of small x evolution, based on the Balitsky-Kovchegov equa- 
tion [39], W[p] can also be modelled by a Gaussian [37], albeit non-local, with fi'^ — > /i^ (^x)- 
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topology corresponding to difFractive scattering off quarks localized in a region 
1/Qs in the nuclei. There is a rapidity gap between the produced particles and 
the scattered quarks in one of the nuclei or the other. The quarks on opposite 
sides of the cut may be localized at different transverse positions^. The other 
two terms with leading contributions have the structure of an interference graph 
depicted in fig. 7(b). 

Of the four remaining terms, two are suppressed respectively by additional 
powers of p and q and two give i5-function contributions for p = ±q. The delta 
function terms are also suppressed relative to the terms we keep at large p and 
q, and would give a contribution not localized in the transverse coordinate, so 
that they would give a flat background once flow effects were included. Both 
of these types of contributions are represented in fig. 8 and are computed in 
Appendix B. We also address there subtleties related to the possible singular 
nature of these contributions. 

The four leading contributions are computed in Appendix A and are shown 
to be identical. So, multiplying eq. 31 by a factor of four, we obtain the leading 
two particle Glasma correlation to be 



S^ (gVj« ^JV,^(A^,^ - 1) 

(2^)6 g4Q2 p4^^4 



•^(P' 9) = 7KZ^6 —TTri Za-a • (17) 



The relation of g^ [i^ to Qs can be quantified numerically by computing Wilson 
line correlators in the nuclear wavefunction. A careful comparison [45] (see also 
Ref. [46]) gives Qs w 0.57 g^/j,^. It is instructive to express the result in eq. (17) 
in terms of the inclusive single gluon spectrum. This is the Gunion-Bertsch [41] 
result and has been been computed previously in the CGC framework [4,11,42] 
to have the form 

dN \_S^ {9^^l,Y N,{N^ 1) ^^ /p^\ ^^g^ 



.dyp(Pp^/ 87r4 g2 ^4^ 

Substituting eqs. (6), (17) and (18) in eq. (3), 

n / dN \ / dN \ 

we find k ^ 4. Identifying the theoretical error on this number is difficult at this 
stage; it requires a numerical computation two particle correlations by solving 
classical Yang-Mills equations as was previously performed for single inclusive 
gluon production [5-7]. 

As we discussed previously, what is measured experimentally is the ratio 
of correlated pairs to the square root of the product of mixed pairs, defined 
as Ap/^pref , where Ap is the difference between the density of measured pairs 
minus mixed pairs and prof denotes the product of the density of mixed pairs [43] . 



^This only makes sense [40] within the framework of an effective theory where one is not 
sensitive to diffractive excitations over some typical transverse scale of size 1/Qs- 



12 




-tr 



eOMAMMMlKl(l(L(L<LI!.(LIL 



P 



4 i 

ii 




^11 

j..C!.l!.«.(!.(L(L0.(BMaa5.(l5.Cl(llLll^ | 



=1= 



^ 



Figure 7: Topology of color correlations. These contributions are detailed in 
Appendix A. The upper (lower) contractions are pairwise contractions of pi 
{p2)- a) Top Figure: Single diffractive contribution to the classical two particle 
correlation. Two gluons are emitted from the same quark line in the amplitude 
and likewise in the complex conjugate amplitude. This diffractive emission how- 
ever occurs at different spatial positions for the sources, which are localized in 
a transverse area of size 1/Qs- See text. There is an identical contribution with 
pi ■*^ P2- b) Bottom Figure: Interference contribution where the transverse 
positions of the interacting quarks are switched in the complex conjugate am- 
plitude for p2 while they are the same for pi . There is an identical contribution 
for pi -^ p2- 
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Figure 8: Topology of color correlations. These contributions are detailed in 
Appendix B. The upper (lower) contractions are pairwise contractions of pi 
(p2)- a) Top Figure: This graph is similar to the single diffractive contribution 
in fig. 7 (a), except, the emission of the particle with momentum p in p2 is 
from different quarks in the amplitude and the complex conjugate amplitude. 
Likewise for the particle with momentum q, where the order of quark lines is 
reversed. There is a similar contribution for pi <-> /92- b) Middle Figure: Double 
diffractive contribution. Two gluons are exchanged from a single quark line in 
both pi and p2- As previously for the single diffractive contributions, the quark 
lines are at different transverse positions in the amplitude and the complex 
conjugate amplitude, c) Bottom Figure: Non diffractive contribution where all 
the quarks are swapped in the complex conjugate amplitude. 
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In our language, the corresponding quantity in the Glasnia can be expressed 
simply as 

/ dN_\ 

^P ^ ciry a) ^^ - ^" (20) 

\ dyp pdp d(t>p I \ dyq qdq d0, / 

where (dN/dy) was defined in eq. (4) and K^^ — kk' k, 4/13.5 ^ 0.3. Note 
however that our computation was performed for large p_L,q±_ 3> Qs while we 
are interested in the p±,q± < Qs region. While we expect the structure of 
eq. (20) to be quite general, as mentioned earlier, we cannot trust the accuracy 
of this prefactor. We will therefore only assume it is a number of order unity to 
be determined by a more accurate numerical computation. 

The expression in eq. (20) is very interesting because it is independent both of 
the rapidities^^ yp and yq of the particle pairs as well as of their azimuthal angles 
(f)p and (f>q respectively. It confirms our picture of flux tubes of transverse size 
1/Qs stretching between the two nuclei (as shown in fig. 4) emitting particles 
isotropically, with equal probability, along their length. This is not the full 
picture though. In the high parton density environment created in central heavy 
ion collisions, the pressure created by interactions among those particles leads 
to collective radial flow. The particles emitted by the Glasma tubes will also 
experience this collective flow. As we shall now discuss, this coUimates the 
relative azimuthal distribution of the pairs. 

We begin by introducing the rapidities of the particles in the direction of 
radial flow (the particles azimuthal angles 4>p^q are defined with respect to the 
radius of the point of emission), (p^q = — ln(tan(0p^g/2)). (It is important to 
note that the two particles will experience the same radial boost only because 
they are localized within l/Qs of each other in the transverse plane-and there- 
fore lie within the same fluid cell.) Expressing the angular distribution (which 
is independent of (/)p and (j)q) in terms of these variables, and boosting it in the 
direction of radial flow, one obtains^^ 

^(P'9) °^ U7~\ T77~^ ^^* U7 T~\ T77 TT ■ (^1) 

cosh(Cp) cosh(Cg) cosh(Cp - Cb ) cosh(Cg - Cb ) 

Here, C^ is the rapidity of the boost and is given by tanh^^ = Vr, where Vr 
is the radial flow velocity. Defining $ = ((j)p + (t)q)/2 and Acf) — (pp — (j)q, and 
re-expressing the boosted distribution in terms of <I> and A(f>, one can re-write 
eq. (20) as 

Ap ,, , , X K„ cosh Co cosh Co ,„„x 

-^{<^,Acf,,yp,yq)^^-^— -^^ -^ — . 22 



^"Quantum corrections, not considered here, will introduce a modest dependence on rapidity 
over scales Ay ^ aj . 

^^The hyperbolic cosines in the denominator come from the Jacobian of the change of 
variables. 
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Substituting cosh^p = 1/ sin (f>p and sinhCp = cos (f>p / sin (j)p, we finally obtain 

J VProf 



K., f d$ 



asiQs) J [coshCB-cos($ + ^)sinliCB] [coshCB-cos($-3^)sinhCB] 

— 71 

K^ 27rcoshCB 



asiQs) cosh^ Cb - sinh^ (^ cos2 ^ ' 
In the particular cases of A(/) = or A0 = tt, the result is 



(23) 




as(Qs) 



27r7^ 

(24) 



where 7^ = cosh^^ is the 7-factor of the boost. Hence, the amplitude of the 
peak, relative to the pedestal, is given by 

A = K,^-l^^^. (25) 

The factor K^ equals 2ttKj^ times the fraction of detected mini- jets "'^2. 

From blast wave fits to the RHIC data, the PHENIX collaboration [44] has 
extracted the average transverse velocity (Vr) as a function of the number of 
participants in a heavy ion collision. To estimate the centrality dependence of 
the coupling^^, as{Qs), we note that the square of the saturation momentum 
<5s — 1-1.3 GeV^ for central Au-I-Au collisions at full RHIC energy, decreasing 
like iVpart [49] towards peripheral collisions'^. The magnitude of A fixes Kn « 
0.6, or Kj^ ^ 0.1, in the ballpark of our simple earlier estimate'^. The resulting 
^(TVpart) is compared to preliminary STAR data in fig. 9. 

The angular width of the correlation function is not reproduced very well 
by the simple "radial boost" model; the integral from eq. (23), without any 
prefactors, is shown as a function of A0 in fig. 10. The width narrows to ~1 
radians only for boost rapidities Cs ~ 2, corresponding to boost velocities of 
about 0.96. To improve the agreement with the measured angular distributions 
one should probably account also for absorption of high-p_L particles by the 
medium [28]. 



^■^Mini-jets from the center of the nuclei which experience little transverse flow or mini-jets 
emitted close to the surface are not detected in the ridge. An accurate estimate of the detected 
mini-jets relative to the total number requires detailed models of the nuclear geometry and 
flow profiles. 

^•^We determine the running coupling from the one- loop QCD /3- function with /3o = llA^c ~ 
27V/ = 27, assuming Aqcd ^ 200 MeV. 

^■'Thc dependence of Q^ on centrality is in fact more complex (we refer to refs. [50]) but 
the simplified form Q^ ~ ^pLit '^ sufficient for our present purposes. 

^^Note, as previously mentioned that only a fraction of the jets are detected, so this number 
is closer to our result that one would naively anticipate. 
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Figure 9: Evolution of the amplitude of the "ridge" (peak at A^ = relative to 
pedestal at A0 = tt) with the number of participants. The preliminary STAR 
data shown is from ref. [43] . 
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3 Discussion 

In the previous section we outlined a novel classical contribution to long range 
rapidity correlations in the Glasma. Boost invariant Glasnia flux tubes of size 
1/Qs in the transverse plane of the nuclei produce particles isotropically with 
equal probability along the length of the flux tube. The common transverse ex- 
pansion experienced by particles in a flux tube coUimates the particles forming a 
ridge like structure in the process. Our result provides a qualitative explanation 
of several features of the ridge [23,24,43,35] in addition to those we discussed 
above. We outline these below. 

• The classical two gluon contribution identified here is qualitatively differ- 
ent from the usual jet mechanism of two particle correlations. It has the 
same structure as single gluon production. This strongly suggests that 
the particles forming the ridge have the same particle composition as the 
bulk spectra. 

• The large classical two gluon contribution qualitatively explains the large 
multiplicity of particles in the ridge relative to particles in the bulk; in- 
deed, we obtain a ratio of order unity naively but this needs to be corrected 
for the fraction of "mini-jets" that are not detected. Estimating this accu- 
rately requires detailed modeling. However, naively comparing our model 
computation to the experimental result, we can guess that this factor is 
about a third of all mini-jets. This contribution increases with centrality 
and energy because Qs is a function of both. 

• The ridge appears independent of the trigger p±^ and is seen for associated 
particles both with a minimal p±^ cut and for p_L,assoc. > 2 GcV. Our 
result was derived at large p_L , but has a geometrical form that seems very 
general and plausible for low p±^ particles as well. The form of our high 
p±^ , q±^ result suggests that the ridge yield is insensitive to the trigger p_L . 
At high p^, it is unlikely the coUimation in azimuthal angle is provided 
by transverse flow; it may instead be an opacity effect resulting from 
the strong dependence of the energy loss of partons on the path length 
traversed in the medium. 

• The ridge amplitude in our model is determined by Q^ and the average 
transverse velocity (K-)- Both of these quantities grow with the central- 
ity (number of participants) and the collision energy. As the transverse 
particle density also grows with centrality and collision energy, it is plau- 
sible that one obtains an approximate scaling of the ridge amplitude as 
a function of the transverse particle density rather than with the energy 
and centrality separately. 

• A study of unlike and like sign charge pairs by STAR [48] shows that for 
small A?7, the signal for unlike sign charge pairs dominates over that for 
like sign pairs, as one would anticipate from jet fragmentation. However, 
for large Ary, the signal for like and unlike sign pairs is the same. This 
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is consistent with a picture of correlated charge neutral sources, such as 
gluons in our case, widely separated in rapidity, emitting opposite charge 
pairs. With such emissions, like sign pairs and unlike sign pairs are equally 
likely to be correlated, just as seen in the data. 

A more quantitative comparison to the ridge data requires more detailed 
numerical computations. Nevertheless, it is encouraging that qualitative fea- 
tures of the ridge are consistent with the formation and flow of Glasma flux 
tubes. We should also note that our picture has features in common with those 
deduced [47] from an analysis of the azimuthal correlation data at RHIC. 
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Appendix A 

We shall first compute the "single diffractive" diagram in fig. 7(a). Performing 
the contractions among the p's corresponding to this graph, we can write this 
contribution to eq. (12) as 

.^(1) ^ (2^f il\ S'^'^i'^^'S^' 5{kij_ - k2±)S{k3J_ - k4±) 

X(5(q_L + P± - fc3_L - ki±)5{qj_ + p± - fc4_L - fe2±) , (26) 

where the superscript denotes the contribution of this particular graph. Substi- 
tuting this expression in eq. (10), we obtain correspondingly, 

^'''(^'^)^^4^)^^^^^/'''^^^^'''(^^''=^^)-^*''('^^'-'=^^)' 

(27) 
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where S± is the transverse overlap area of the nuclei and^® where 

Studying the structure of eq. (29), it is clear that the terms with the fewest 
powers of p_|_ , q_|_ in the denominator (and therefore the largest contribution to 
C} ) are those with the fewest powers of fcu^ in the numerator. One then finds 
simply that 

^W(p^,fci^)^«(q^,-fci^)^^^-^. (30) 

Note that the angular dependence on the relative angles of feij^ with pj^ and 
q_i_ cancels between the leading longitudinal and transverse contributions to 
eq. (29). Substituting this result back into eq. (27), this gives 



4(2^)6 g4g2 p4^qA_ 






where we have used the infrared cut-off fcmin — Qs, which is the saturation scale 
scale signifying the onset of non-linear contributions that soften the infrared 
gluon distributions in the CGC. 

We now turn to the "interference diagram" diagram in fig. 7(b). Again, 
performing the contractions among the p's corresponding to this graph, we can 
write this contribution to eq. (12) as 

.^(2) = (2^)« Ai« ^b''^^-^- S{k^^ - k21_)S{k^^ - fc4±) 

^S{pj_ - q_L - fc3± + fc2±)<5(q_L -P± + fci-L - fe4±) , (32) 

where the superscript denotes the contribution of this particular graph. Again, 
plugging this in eq. (10), we obtain, 

(33) 
The leading contribution to the integrand is the same as in eq. (30), and we 
obtain 

a"'(p,.)=J-, '<%''-"'T'r^ (34) 

4 {2TTf g4 Qi p\ q\ 

which is identical to eq. (31). 



^''In general, 

-4(Px , fcix , fc2x ) = To—, r — S2.2 , r — ^ ' (28) 

HAPI- - fcix) fc2x(P± - k2l_Y 
where the product of Lipatov vertices is given in eq. (11). In eq. (27), we use the notation 

^(Px,fcix,fcix) = ^(i)(px,fcix)- 
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The other two leading "single diffractive" and "interference" contributions, 
which we label by C^^-* and C^^-* respectively, have a structure such that 

C(3) < — > C(i) for pi < — > p2 

C(s' < — > C(2) for pi^^ P2, (35) 

so it is not surprising that these terms give contributions that are identical to 
eqs. (31) and (34). 

Appendix B 

We now turn to the topologies discussed in fig. 8. Consider first the topology 
shown in fig. 8(a). We obtain in this case the contribution to eq. (12) as 

.^(5) = (2^)8 /iS fbgbb^c-c^cc ^(^^^ ^ ka^)5{k2± + fe4±) 

X'5(P_L - 9± - fc2± + k3±)S{q_^ -P±- k4± + ki±) , (36) 

where the superscript denotes the contribution of this graph. Again, substitut- 
ing this expression in cq. (10), we obtain 

^ ^^'^^-64(2^ 2 ^^^^ 

X / (fkj_Aipj_,kj_,pj_ -q_L - kj_) A{qj_, -kj_, q - p + kj_) , 

where, from eq. (28), 

4 

A{pj_,kj_,pj_ -q_L -fe±) = T 

PI 

(P_L - ^-l) • kj_ipj_ -qj_- k±) ■ (q_L + kj_) 
+ (fc_L X (p_L - fe-L)) • ((P± - g± - fe-L) X (q_L + '^-l)) 



(37) 



/ci (p^ - fe^)2 (p^ - q^ - fci)2 (qr^ + fc^)2 



,(38) 



and likewise, a similar expression for A{qj_,—k±,qj_ — pj^ + k±). Unlike the 
expression in eq. (30), this expression docs not have a simple expression in the 
infrared. It is however suppressed by additional powers of pj^ and Qj^ relative to 
the leading terms, allowing us to ignore this sub-dominant contribution. There 
is another contribution, that we denote by C'-^^ obtained by exchanging pi ^^ p2 
in fig. 8 (a) that is similarly suppressed. 

Finally, there are the diagrams discussed in fig. 8 (b) and (c). The former 
has the form of a "double diffractive" contribution. We denote its contribution 
to eq. (12) by 

.^(4) ^ {27Tf m' S^'S^'S'^S^' Siki^ + fe3±)'J(fc2± + fe4±) 

xd{pj_ + g_L - f^2± - k4±)d{qj_ +p_L - fci-L - fc3±) • (39) 
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From the structure of this equation, we will have 

X / rf2fci_L d2fc2± ^(p_L, fel-L, fe2±) ^(-P±, -fel±, -fc2±) , (40) 

where the form of yl is obtained from eq. (28). One obtains a similar expression 
for fig. 8(c), where instead of a (5- function in pj^ + q^^, one gets instead a 6- 
function in pj^ — q_|_ . These (5- functions will be smeared out by re-scattering and 
are suppressed relative to the leading terms in Appendix A. 
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